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ABSTRACT 

We derive exact expressions for the degree of lineal polarization over a resolved or integrated stellar 
disc due to resonance scattering and the Hanle effect from a dipolar or quadrupolar distribution of 
magnetic fields. We apply the theory of scattering polarization within the formalism of the spherical 
tensors representation for the density matrix and radiation field. The distribution of linear polarization 
over the stellar disk for different configurations of the magnetic field is studied and its topology 
discussed. For an unresolved dipole, the resulting polarization can be expressed in terms of just three 
functions (of the inclination angle and effective dipole strength), that are calculated numerically and 
their behaviour discussed. Dipolar and (aligned) quadrupoles are considered in some detail, but the 
techniques here — in particular, the extensive use of the spherical tensor formalism for polarization — , 
can easily be applied to more general field configurations. 

Subject headings: line: formation — polarization — scattering — stars: magnetic fields 
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The scattering of light in a plasma produces linear polarization in spectral lines and continuum. As a consequence, 
I ■ the spectrum from a stellar atmosphere is linearly polarized towards the stellar limb. A remarkable example of this is 
O ' the linearly polarized component of the Fraunhofcr spectrum observed close to the solar limb, also called second solar 
spectrum for its radically different nature and str ucture with respect to the intensity spectrum (jStenflo et al.l fl983: 
c/3 . iStenflo fe Kellerlll997t[Gandorferi[2000ll2002l [20051) . In other (unresolved) stars, its observation is difficult because the 
scattering polarization spectrum cancels when average dd over the whole stellar disk, its observation requiring imperfe ct 
cancellations due to, for example, eclipses or transits (jLandi Degl'Innocenti et aT1ll988l : iCarciofi fc Magalhacs 2005). 
The presence of a weak magnetic field partially removes the degeneracy of the magnetic sublevels of the scatterers, 
■ which perturbs the scattering process and hence, the polarization pattern (Hanle effect) ()G. Moruzzi fc F. Strumial 
f*"*- ! 119911 ). By weak we mean that 0.1 < vj/r < 10, where the vl — 1.3996 x 10 6 £? is the Larmor frequency (expressed 
OO 1 in s" 1 ) corresponding to the magnetic field B (in G), and r the characteristic time for the scattering process - 
' for resonance scattering with unpolarized lower level, t ~ l/A u e, A u g being the Einstein coefficient for spontaneous 
. emission in the transition^. Typically, this is a regime of fields weaker than the characteristic regime for the Zeeman 
Q\ • effect (vl ~it Afli/; Arjv being the Doppler width of the spectral line). The possibility of exploiting the Hanle effect for 
\ diagnosing solar magnetic fields in magnetic field regimes and topological configurations not easily accessible through 
CNJ . classical Zeeman techniques is one of the reasons why the second solar spectrum has been so activel y researched in 
t— ( ■ recent years (see lStenfloNl982l iCasini fc Landi DeglTnnocentil 120081 [TVu iillo Bucno 2009: lStenflo1[200l and references 
' therei n) . Moreover, the possib i lity has arisen that the Hanle effect could also be applied to diagnose stellar magnetic 
. Jh '■ fields (jlgnace et al.lll995l Il997t lNordsieck|[200l1: Ilgnacell200ll and references therein). 

The presence of a global magnetic field over the stellar surface brakes the symmetry limitations of the pure scattering. 
Hanle effect signals depend on the orientation of the magnetic field with respect to the line of sight (LOS) and on 
, the geometry of the scattering event (i.e., the relative orientation of the magnetic field with respect to the stellar 
surface). It is then possible that scattering polarization signals appear even at disk center (for inclined magnetic fields 
there) and not only towards the stellar limb, and perfect cancellation when averaging over the stellar disk is extremely 
unlikely (but for very special field configurations). Here, we derive exact expressions for the intensity and polarization 
produced by scattering polarization and th e Hanle effect in the presence of a poloidal distribution of fields (see also 
lLopez Ariste et aD 1201 It llgnace et al.ll2011| ). We will consider with some detail relatively simple distributions of the 
magnetic field (a dipole, an aligned quadrupole) over the stellar disk, though more general configurations could be 
studied along the same lines. Stellar dipoles have been used to study the polarizati on due to the longitu dinal Zeeman 
effect from stars fe.g„ Schwarzschildl Il950l ). th e transversal Zeeman effect (e.g.. lLandolfi et al.lll993l ). synchrotron 
radiation from stars ( Thornd 19631) and planets (|Chang fc Davislll962l ). 



The next section summarizes theory of the Hanle effect in a resonance transition, within the spherical tensor com- 
ponents formalism of the density matrix. The theory is then applied (Section 3) to calculate the linear polarization 
emission in a stellar disk with a poloidal distribution of fields, for which extensive use of the spherical tensors formalism 
is done. The distribution of polarization on a resolved disk is discussed in Section 4; the polarized emission by an 
unresolved oblique rotator in Section 5. The discussion is extended to quadrupolar fields in Section 6. 

This is just the slowest timescale in the absorption-reemission process. The faster absorption time scale T a j, s = X/B( U J, with Bg u 
the Einstein coefficient for absorption and J the mean intensity, is irrelevant here since the lower level will be assumed to be unpolarized 
(Section 2), thus playing no role in the polarization of the scattered light. 
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Fig. 1. — (a) We consider a Cartesian reference system with the z axis along the line-of-sight (LOS), and the stellar dipole lying in the x-z 
plane, at an angle i with respect to the LOS. The reference polarization direction ei (positive direction for Stokes-Q) forms an arbitrary 
angle \ with the x axis; correspondingly, e2 is at an angle \ with respect to the y axis, (b) In the stellar reference system, z' lies along the 
dipole, and y' = y. At an arbitrary point on the stellar surface with longitude I from the x'-z' plane, and colatitude A, the magnetic field B 
lies in the meridian plane (poloidal field), and is inclined an angle 8g with respect to the vertical, (c) Averaging over the projected stellar 



-dxdy = Jq n pdpda (d) In an oblique rotator, the magnetic 



disk can be done in a number of different ways: ff d£ = f^} R J^^— j- 

axis e is inclined at an angle \E" with respect to the rotational axis oj which, in turn, is at an angle j to the LOS. We choose the reference 
direction ei for polarization, along the projected rotation axis (i.e., in the x"-z plane formed by u> and the LOS). 

2. SCATTERING LINE POLARIZATION AND HANLE EFFECT 

Neglecting the Zeeman splitting between a and tt components of the line profile (Hanle effect regime), the emissivity 
in the Stokes parameters /, Q, an d U along a given line-of-sight (LOS; see Figure fa) in a transition with a polarized 
upper level can be expressed as (jLandi Pegl'Innocenti et all If 9901 : lLandi Degl'Innocenti fc Landolfil I2004L hereafter 
LL04): 
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where e Q — j^A u tN \J2J U + l(j> v (A u i is the Einstein coefficient for spontaneous emission, N the number density of 
atoms, and cp v the absorption profile), Wj^j e is a numerical factor depending on the total angular momentum J of the 
levels involved in the transition (see lLandi Degl'Innocenti fc Landolfil (|2004f ) for explicit values), and pq [K = 0, 2J; 
Q = —K, ...,K) are the spherical components of the density matrix (|Fanoi ri957) of the upper level of the transition 
(defined here taking the quantization axis along the LOS), which can be expressed explicitly in terms of its real (pq) 



and imaginary (pq) parts: pq = pq +ipq. (In particular, p^ is f j \JU U + 1 times the relative population of the upper 

level; p\ a measure of the imbalance of populations between the magnetic sublevels with a well defined projection Mj 
of the angular momentum on the LOS). 

In the magnetic field reference system (quantization axis along the magnetic field), the density matrix satisfies 
(jLandi Degl'Innocenti et~allll 990L LL04): 



[l + *W(l- e )] 



2hv 3 2J t + 1 
c 2 



9% 



1 



f + irg 



(4) 



where 5^ = D( K ^/A u e, is the elastic collisions depolarizing rate, e = C u i/(A u t + C u i), C u e is the collisional 

de-excitation rate, T — j(l — e)/[f + S K (1 — e)], 7 = 2-Kv^g u j A u i, = f.3996 x 10 6 B is the Larmor frequency 
in s _1 , B the magnetic field in G, g u the Lande factor of the upper level, and B v is the Planck function. The 
radiation field tensor spherical components Jq (u) {K = 0, 1, 2; Q = —K, K) are averages over all directions of the 

radiation field illumination incident upon a point (jLandi Degl'Innocentil[l984l LL04). In particular, Jq {v) = § rl^x 1 



(dSl = sin#d(Mx' = d/zd%' is the element of solid angle), Jq{v) — 



^ )a ndJ » = /f ^ 



[(3/i 2 



Din 



3(/x — l)^^^^'] are the only non- vanishing components if the radiation field is axially-symmetric and independent of x' 
(e.g., in a piane-parallel medium or in a spherically symmetric atmosphere). The radiation field tensors are averaged 
through the absorption profile in Equation (|4]), Jq = J dv(j} v Jq{v), in accord with the complete redistribution 
hypothesis (the frequency of the incident photon is redistributed, collisionally or otherwise, and unrelated to that of 
the incident scattered photon). Hence, Jq = 0, since there is no net circular polarization in the medium. 
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3. DISK DISTRIBUTION AND INTEGRATION OVER THE STELLAR DISK 

A rotation R (Euler angies (a/37)) fr° m an "old" reference system to a "new" one transforms the spherical components 
of the density matrix (the "*" symbol stands for complex conjugate), and radiation field tensors according to the 
following rules 

Q' Q' 

where Vpg^afij) = dpg(/3)e _I ^ Pa+< 5 7 ' ) is the rotation matrix and dpg the reduced rotation matrix (jBrink fc Satchlerl 
EMI). Hence, J K Q = Y Jq [J K q \ ■"d-q-Q^B) in Equation (jH), where [Jl q ] v are the values of the radiation field tensors in 
a reference system with the quantization axis along the local vertical direction, and 9b is the inclination of the magnetic 
field with respect to the local vertical (Fig. lb). Multiplying Equation (gj) by V 2 QQ , (0, -Q B , -£)* and V 2 QIQ „ (0, -i, 0)* 
(0 b = Ob + A, A is the colatitude with respect to the dipole; see Figure lb), and summing over Q and Q' , we get the 
density matrix in the reference system with the quantization axis along the LOS and the dipole vector contained in 
the x-z plane. After some index renaming: 

(i-^iEC-^ra- E T-^<Q'(^)4Q^-e B )4, Q (-z)e-^" 

q Q'Q" H 

+ eB lJ 8 K0 8 Q0 . (5) 

Equation §5§ is in the xyz reference system (see Figure 1), and can thus be directly plugged into Equations (dJ-Q. 

Computing the exact Jq components at all points in a plane-parallel, spherical, or multidimensional atmosphere 
is a challenging numerical problem (e.g.. Ilgnace et aflll999t PTru iillo Bu eno fc Manso SainzlU99 9; Nagendra 2003J. In 
an optically thick atmosphere, Jq ^> Jq 3> |Jq-^oI' which means that the radiation field is only weakly anisotropic 
and that azimuthal symmetry breaking effects (due to horizontal inhomogeneities or magnetic fields) are negligible 
(|Manso Sainzll2002t IManso Sainz fc Truiillo Buendl20"Tl[) . Thus, neglecting Jqm elements in Equation J5]), for pQ-. 

[i + s^(i ^ 2J ^^==pI = (i - ^» E T^AQ'iWl^i-QB^i-i). (6) 

Q'Q 

As a consequence of the weak anisotropy, the polarization in an optically thick medium is low (Q/I.U/I <C 1), 
and the anisotropy of the radiation field is determined by the center-to-limb variation of the intensity J5 2 = 
§ ^ ^7g[(3/x 2 — 1)/ MX '- J q can be estimated from the emergent intensity (e.g., IManso Sainz fc Landi Degl'Innocentil 

20021 IManso Sainz et all [20061. Thus, for example, assuming a linear limb-darkening law for the intensity in the 
spectral line 

I = J(°>(1-u + um), (7) 

wher e \i = cos 9 and it the limb-darkening coefficient, which is known from observations or numerical modeling 
(e.g., iPierce fc Slaughter! U977t lAllende Prieto et al.|[200l iQuirrenbach et al.lll99l iNordgren et al.lll999t Ivan Hammel 
1993: Claret 2000), or otherwise, can be left as a free parameter. For the darkening law in Equation ([7]), [Jq] v = 

JC°3(1 - u/2)/2, and [Jg] v = j(°>u/16v^. 
For a dipole, the distribution of magnetic fields on the stellar surface is 

B = -^[e-3(e.r)r], (8) 

where Bd is the magnetic field strength at the poles, and e is the unit vector directed along the dipole. For a point P 
over the visible disk whose position is indicated by the unit vector r (see Figure lb), Equations (j6|) can be written as 

[1 + *<*>(! - e )]^^^= pi = (1 - eWZMZMpd, T), (9) 

[1 + - e)] 2k f 2 7 ^±J= p 2 ±2 = (1 - e)w^ Je [J^Ug P {i, T) ± ih P (i,T)\, (10) 

where we have introduced the following functions of the dipole inclination i, strength T^, at a given point P over the 
stellar disk: 

f P (i,T d ) = E 1 + r 2 Q/2 rf OQ'(^)4'Q"(- e ^)4"o(-^)[cos(Q^) - rQ'sin(Q"l)], (11) 
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g P (i,T d ) = ]T T -^^rfoQ'(^)4'Q"(- e s)4"2(-*)[cos(Q^)-rQ' S in(Q^)] ! (12) 

Q'Q" + W 

= - £ i 2 rfg Q , (g B )d 2 QIQ „ (-6 B )4„ 2 (-z) [LQ' cos(Q"l) + sin(Q"l)] . (13) 

Q'Q" W 

In Equations (ITT1) - (|13I) . the dependence on P (and also on i), appear implicitly through (see Figure la-c) 



r»=I^[l + 3(e.rn cosA = e .,, C0S 6 B = - /f ^= , cosl= ^ - ^ ^ , 

where $7 is the unit vector directed along the LOS. In the saturation regime (Td — > oo), only the Q' — terms remain 
and 

f P (i, oo) = d 2 00 (6 B ) dl Q '- (-6 B )4„ (-i) cos(Q"^), (14) 
Q" 

ffp(»,cx>) = 4,(fl B )X;^g-(-efl)4-2(- i )coB(Q''^), (15) 

Q" 

h P (i, oo) = -dg (d fl ) £ 4q-> (-©b)4"2(-«) sin(Q^). (16) 

Q" 

Equations (|9j)- (fl"0)) are expressed in a common reference system — with the quantization axis along the LOS — , and 
they can be averaged over the visible disk 

[l + fW(l - e) ]^_^+lg = (i- e) W % t [J*] v f{i,T), (17) 

[i + ^(i - e)]^^±J= P 2 ±2 = (i - e )«,w ji [^].&(i,r) ±ifc(i,r)], (is) 

where we have introduced the following functions of the dipole inclination i and strength F^: 

f(i,T d ) = i // dS E 1 + L o A ( 9 B)dl Q ^-eB)d^ (-i)co S (Q''e), (19) 

* JJ Q'Q" ^ 

g(hTd) = ^2 // ds E 1 + r2 Q ,2 rf 0Q'(^)4'Q"(- e ^)4"2(-»)cos(Q^), (20) 

= ^//dS E ! ~r?Q/2 ^Q'( g g) d Q'Q"(- e ^) d Q"2(-0cos(Q^). (21) 

In Equations (fT9 |) - ([2~T|) . the integral extends over the projected disk (see Figure lc) and i?* is the stellar radius. In 
passing from Equations (ITT1) - (|13I) to Equations (|19l) - ([2"Tj) . we have taken into account that two points at the same 
colatitude A, and longitudes ±£, have the same values of magnetic field strength (r) and inclination (9b) — hence, the 
same <d B too. Therefore, averages with sm(Q" t) cancel out. 

The functions f(i,Yd), g(i,Yd), and h(i,Yd) are evaluated numerically and are shown in Figure 2. By construction 
they do not depend on the azimuth of the dipole with respect to the line of sight. From this fact and the properties 
of the d 2 QQ, rotation matrices in Equations (ITT1) - (jl3l) . the the following symmetries follow: f(i,Td) = /(tt — i,Yd), 

g(i, Yd) = g(^ — i, Yd), h(i, Yd) = —h(ir — i, Yd). It is also clear from symmetry arguments that in the absence of fields, 
the integrated linear polarization vanishes — hence, g(i,0) = h(i,0) = — , while /(i,0) = 7r/43- In the saturation 
regime (Td — > oo), h(i, oo) = 0, and 

/(*, oo) = i J J dS d 2 00 (9 B ) d 2 0Q < (-e B )4„ (-i) cos(Q"£), (22) 

Q 

9(i,°°) = jpjj dS^ (^)E d 0Q"(-©B)4" 2 (-*)cos(Q"£). (23) 

Q' 

Consequently, in saturation, p\ = p\ = 0, and ejj = 0. Light is thus polarized along the dipole and the maximum 
amount of linear polarization will be obtained at an inclination i = 90° . 

2 Clearly, f(i,0) = /(0,0). Now, making use of d 2 QIIQ {i = 0) = Sq/iq, and T, N d J MN {a)d J Np {l3) = d J Mp (a + j3) (e.g., IBrink fc Satchlerl 
[1963 . LL04), then /(0, 0) = //<Gdg (-A) = tt/4. 
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Fig. 2. — Functions f(i,T d ), g(i, F d ), and h(i,F d ), defined in Eqs. l|f 9jl - l|2f jl . for different values of the inclination angle i of the dipole 
to the LOS. For i > 90°, f(i,T d ) =_/(180° - i,T d ), g(i,T d ) = §(180° - i,T d ), and h(i,T d ) = -fe(180° - i,T d ) Note that f(0,T d ) = 
/(90°, Tj) = g(0, F d ) = 0), and that h(i, F d ) reaches its maximum at i £S 52°. 

The emergent Stokes parameters along the LOS (z-axis) from a semi-infinite atmosphere^ are 



/>00 />OG />OG 

I = Si exp(— r)dr, Q = Sq exp(— r)dr, U = I Su exp(— r)di 
Jo Jo Jo 



(24) 



where 5"/ = e//ry/, S Q = e Q /r]i, eu/r)i, Vi = 7 7/ inc + Vi° nt ' 7 ?/ n ° = ^B tu N t <j> u (Ben is the Einstein coefficient for 
absorption), and dr = —rjidz is the element of optical depth, which is measured from the stellar surface (r = 0). The 



3 The atmosphere cannot be considered semi-infinite very close to the extreme stellar limb. We neglect those effects here. 
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integrals are approximated by the value of corresponding source functions Si at r = 1, which we evaluate from the 
emergent center-to-limb variation. The emergent fractional polarization then reads 

9. = £0 = e Q n ° aS ^ e X L = e JL = e 'u C aS 'i ne (25) 

I ei e l } nc aSf ne + B v ' I ej e 1 / 110 aS*' inc + ' 1 j 

where a = rj\ lnc / 7]j° nt , and £?„ = € c ° nt / rff ont and S 1 } 1110 = e l j nc / ri\ lnc , are the source functions for the continuum (Planck 
function) and line, respectively. At the line core of a strong line, a ^> 1 and Q/I = fg 10 / 6 /™! U/I — ey ne /ej me ; for a 

weak line, Q/J = e^/e^a^^l - a^T + ■•■), = ^V^a^C 1 " «4r + -)■ 

If the emissivities e, (see Equations HJ-©) are evaluated from Equation ([S]), then Equations (|2"5")l express the 
fractional polarization at a given point on the resolved stellar disk; if they are evaluated from Equations (|9|)- ([T3)) . 
then we get the polarized fluxes Ti, J-q, and J r u, integrated over the stellar disk. For symmetry reasons, e C Q° nt does 
not contribute to the integrated fractional polarization in Equation (|2"5")) . though Rayleigh and Thompson scattering 
polarization should be included if the stellar disk were resolved. If the reference direction for positive-Q is taken along 
the dipole projected on the plane of the sky (i.e., \ — in Equations ([2])-([3])), then Ti, J t q, and Txj are univocally 
determined by the functions /, g, and h, respectively. In particular, the symmetries discussed above for /, g, and h 
apply to J- 1, J-q, and J-jj ■ Finally, note that the not only the linear polarization, but also the total flux J-j, is modulated 
by the presence of a global dipolar field on the stellar surface (through f(i,Td))', this modulation is maximum for a 
dipole aligned with the LOS and minimum for a transversal dipole. 

We note in passing, that the maximum amount of fractional polarization is obtained, in the absence of magnetic 
fields, at the stellar limb; from Equations (H])-© 



,na X = V^+TP J 3iy 2 [J Q 2 ], 

PL I max 2V2[l + 5(l-e)][JS + T ^B,}-W 2 [J§} v ' [ ° J 

where W<i = {w^jj 2 is the line polarizability. If there is no limb- darkening or brightening, u = 0, [Jq] v — 0, and hence, 
= 0; if the radiation field at the surface is highly collimated ([Jq] v = Jq/V2), and the atmosphere collisionless 
(e = 5^ = 0), the ideal limit of a 9 0° single-scattering event p™ ax = 3W2/(4 — W2) is recovered (cf., Equation (10.26) 
in lLandi Degl Innocenti fc Landolfil 120041 ). 

4. RESOLVED DIPOLAR DISK 

It is illustrative to study the distribution of intensity and polarization over the resolved disk of the dipole. This 
problem may also be of practical interest for example, to study the emission of planets at long wavelengths where their 
intrinsic emission dominates over the irradiation of the Sun. 

Using Equations (|9|)-(fT3|) into Equations ([I])-© and then, into Equations ([25]). we obtain the emergent fractional 
polarization at every point over the stellar disk (Figure [3]). We consider a strong line a = 10 3 , with NLTE parameter 
e = 10~ 4 . In the absence of magnetic fields (upper most row in Figure [3]), the intensity shows the limb darkening law 
(Equation ([7]); second column of Figure the polarization increases from center (where it vanishes for clear symmetry 
reasons; third column of Figure [3]) to the limb (p maK ), and the polarization plane is always parallel to the stellar limb 
(see fourth column of Figure [3]). 

Now, we consider a dipolar field with = 2, observed pole-on (i.e., inclined at an angle i = 0° to the LOS; second 
row of Figure[3]) . Due to the rotational symmetry, the disk center is still brighter but the limb darkening law is distorted 
with respect to the non-magnetic case due to the modification of the alignment p\ in Equation ([T]) by the presence 
of magnetic fields. The polarization is no longer maximal at the limb, but in a ring within the disk. This is because 
the magnetic field at the equator is pointing away from the LOS, a configuration that maximizes depolarization and 
rotation of the polarization plane characteristic of the Hanle effect. At disk center the magnetic field is vertical, and the 
Hanle effect does not operate (also clear from symmetry reasons). Half way between these to points, the polarization 
reaches a maximum, while the rotation of the polarization angle varies with the relative inclination of the magnetic 
field to the LOS, hence with the distance to disk center. 

When the dipole is inclined i — 30°, 60°, and 90° (perpendicular) to the LOS (rows 3-5 of Figure[3j respectively), the 
magnetic field at disk center becomes inclined to the LOS, and the symmetry arguments above do not apply anymore. 
Thus, the brighter region on the disk is displaced from the center, and limb darkening depends on the distance from 
the disk center as well as on the azimuth. Actually, in some cases, several bright regions may appear. These bright 
spots are not due to abundance gradients (we are considering a homogeneous atmosphere) but just to a modulation of 
the intensity by the magnetic field, and if not properly accounted for, they could distort the interpretation of Doppler 
imaging techniques (see next section). The behavior of polarization is likewise complex, depending on the distance to 
the limb but also on the azimuth. Only in the perpendicular dipole configuration some symmetry is recovered, with 
the distribution of intensity and polarization being symmetric with respect to the equator. 

Finally, it is interesting to observe the behavior of the singularities of the linear polarization field (points where the 
linear polarization vanishes and hence, the linear polarization angle is undefined) with the inclination of the dipole, 
which is most clearly illustrated in the rightmost panels of Figure [3] In the absence of magnetic fields or for a pole-on 
configuration, there is just one singular point for a^: disk center. With an inclination i = 30° two close singularities 
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Fig. 3. — Distribution of intensity and linear polarization over the disk of a dipole. First column on the left shows, as a reference, the 
absolute value of the magnetic field strength \B\ — which varies by a factor 2 between poles and equator. Columns 2 to 4 show the intensity, 

degree of linear polarization Pl = \J Q 2 + U 2 , and polarization angle aj, (tan(2a^) = U/Q), respectively. Uppermost row shows the pure 
scattering, zero field case; rows 2 to 5 correspond to = 2 and inclinations to the LOS i = 0° (pole-on), 30°, 60°, and 90° (transversal), 
respectively. Crosses mark disk centers; diamonds, poles. 

appear slightly off-center, four at 60° (two close ones in the north hemisphere, one in the south, and one close to 
the equator), and one at each hemisphere at 90° inclination. Now, consider a small closed path around a singular 
path and we follow it say, clockwise. For the first two cases (non-magnetic; pole-on dipole), we change twice color 
(i.e., polarization direction), and as we follow such a path we pass from blue to green to yellow. We arbitrarily 
assign to such a singularity an index +1. Proceeding along similar paths around the singularities for the i = 30° 
inclination, we find the same color (polarization direction) change pattern (blue-green- yellow-...), but every color is 
found only once every turn. We assign such singularities an index +1/2. Three of the four singularities appearing 
in the i = 60° inclination case follow this behavior too, and have index +1/2, but the fourth (the one close to the 
equator) changes once in reverse order (yellow-green-blue-...), for which an index -1/2 is assigned to it. Finally, the 
two singularities in the i = 90° case have indexes +1/2. It is then easy to verify that the total amount of these indices 
is conserved (+1 = +1/2 + 1/2 = +1/2 + 1/2 + 1/2 — 1/2). This is closely related to the topology of singular points 
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Fig. 4. — Periodic behavior of the Hanle effect modulated polarized emission from an unresolved oblique rotator with a dipolar distribution 
of fields. From left to right: intensity fluctuation normalized to the intensity at t = 0, linear polarization p^, azimuthal angle a (tan(2a) = 
T\j ITq), and Tq/Ti (solid lines) and T\j /Ti (dotted lines). Upper and lower panels show two different general behaviors. Upper panels: 
i < 90° at all times (blue line: j = 30°, * = 20°; red line: j = 20°, * = 30°). Lower panels: i fluctuates above and below 90° at different 
times (blue line: j = 80°, \P = 40°; red line: j = 40°, it = 80°). Note that the red and blue lines superpose in the two left columns. 

in polarization fields ()Berrv et al.ll2004l ). It would be very interesting to consider the diagnostic value of such points. 
Just counting them or knowing their approximate location could be used to constrain the global topology and strength 
of the magnetic field. 



5. OBLIQUE ROTATORS 

An oblique rotator is a star whose magnetic axis is inclined at an angle to the rotational axis. As a consequence, 
the polarimetric signal varies periodically as the star rotates. 

Let ^ be the angle between the magnetic dipole (e) and rotation (u>) axes, j the inclination of the rotation axis to 
the LOS, and let's choose the positive-Q direction as the plane containing the rotation axis and LOS (see Figure Id), 
then 

cos«(/) = — sin VP cos / sin j + cosVE'cosj, (27) 
sin X (f) smi(f) = sin sin /, (28) 
cos%(/) sini(/) = sin 'J 7 cos / cos j + cos 'J sin j, (29) 

where / (0 < / < 2ir) is the rotational phase angle (/ = when the LOS, u;, and e are coplanar and thus ordered; 
see Figure Id). Using Equations ([27l) - (|29)) into the expressions (fT7|) - (l2"Tj) and ([T])-([3]), we calculate the intensity and 
polarization with rotational phase (Figure 2]) . 

The curves of brightness and polarization fluctuation fall in one of two broad classes according to their inclination of 
the dipole to the LOS being i < 90° at all times during a rotation, or its changing from i < 90° to i > 90°. In the former 
case, the brightness and linear polarization vary once between maximum and minimum values in one period (upper 
panels in Figure 0}; in the latter, there are two maxima and two local minima (lower panels in Figure 2]). Fluctuations 
in total brightness and polarization are symmetric with respect to the middle part of the rotational phase. This is 
because in every rotation, the two configurations of the dipole at ±f (or, equivalently, at / and 2n — /), show the 
same inclination i to the LOS, their only difference being their azimuth, or equivalently, an inessential rotation of the 
re ference system for the po l arization direction . Therefore, Fx{f) = J-i(2n — /), and Pl(I) = Pl{^ — f) ( m contrast 
to lLopez Ariste et al.|[20iH: llgnace et al.ll201lD . 

There are two different behaviors for the azimuthal angle which is illustrated in the third column of Figure (0]) . 

Figure [5] shows polarization diagrams [Tul^i versus Tq/Ti) along the rotation cycle, for different values of the 
j and 'J 7 angles. These diagrams are computed choosing the reference direction for positive-Q along the projected 
rotation axis, and considering the only Hanle effect polarization. In actual observations, the orientation of the rotation 
axis on the plane of the sky forms an unknown angle ip with respect to the arbitrary reference direction chosen for 
polarization. On the other hand, interstellar polarization is often not negligible, adding an arbitrary amount of linear 
polarization (constant along the line profile) . The former effect amounts to a rotation by an angle 2ip of the diagrams 
in the Q-U space; the latter to a translation. Therefore, the orientation and position of the diagrams may change, 
but their shape remains unchanged, being an intrinsic characteristic of the dipole. Diagrams similar to these have 
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Fig. 5. — Polarization diagrams Tjj /Ti versus J-q/J-j for a dipolar magnetic field configuration with = 2 (solid lines) and in saturation 
(r^ = oo; dotted lines), at an angle to the rotation axis, which in turn forms an angle j with the LOS. Each diagram corresponds to a full 
star rotation. The diagrams show polarization varying between ±1% (grid with steps of 0.5%). We consider u = 0.2, e = 10 -4 , 8^' = 0. 

been obtained for different mechanisms capable of generating l inear polarization in spectral lines, like for example th e 
differential saturation or magnetic intensification mechanism ()Landi Degrinnocenti et aLl Il981t lLandolfi et al.lfl993l) . 
although for more intense magnetic fields. 

The degree of integrated linear polarization seen at ±f (or / and 2n — /) is the same, but the diagrams are 
not necessarily symmetric with respect to the J-jj = lane in general. In the saturation regime, however, the 
polarization signal depends on the direction of the magnetic field, but not on its orientation; hence, J : Q{f) = Fq{— f), 
Fu{f) = —Fu{—f), and the diagrams are symmetric with respect to Txj = 0. 

6. STELLAR QUADRUPOLES 

We consider a distortion of the dipolar field in the form of an additional multipolar component of the global field. 
For the sake of simplicity, we consider a quadrupole aligned with the dipole e. The contribution of the quadrupole to 
the distribution of magnetic fields on the stellar surface is 

B q = -§42(e • r )e + [1 - 5(e • r)>}, (30) 

where B q is the magnetic field strength at the poles. The total magnetic field on the stellar surface is now B = Bd + B q , 
where Bd is given by Equation ([8]). The global field is still a poloidal field, and Equation © applies. From it we may 
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Fig. 6. — Same as Figure 5 but for a purely quadrupolar distribution of fields (q = oo, T q = 2; dotted lines), and a dipolc+quadrupolc 
(q = 1, = 2; solid lines) distribution of magnetic fields. 

derive expressions for the resolved and integrated disk emission as in Equations (|9|)- (fTU| and (fT7| -(fT8 |) . respectively, 
but introducing more general functions fp(i, Td, T q ), gp(i,Td,T q ), hp(i,Td,T q ), and the corresponding f(i,Td,T q ), 
g(i, Td, h(i, Td, where T q is given from B q as after Equation (|4]). The expressions of these functions are formally 
identical to their dipolar counterparts (Equations (fTTj) - (lT3|) and Equations (fT9|) -([2T |) '). but now their dependence on 
the quadrupolar component is implicitly given by the normalized strength and inclination of the local magnetic held 

r 2 = 3[i + 3 (e • rf] + ^[1 - 2(e • rf + 5(e • r) 4 ] + 2T d T q {e ■ r) 3 , 

, = 2(e-r)- g [l-3(eT) 2 ]] 

3 [(1 + q 2 ) + (3 - 2(7 2 )(e • r) 2 + 8q{e ■ r) 3 + 5q 2 (e ■ r) 4 ]V2 ' 

where q = B q /Bd. 

Figure [B] shows polarization diagrams for a dipole+quadrupole configuration (r^ = 2, q = 1), and a pure quadrupole 
(Td = 0, T q = 2). A purely quadrupolar distribution of fields tends to distort the polarization diagrams in an opposite 
sense to the dipolar distribution considered in Figure [5J and consequently, the dipolar+quadrupolar distribution tends 
to be more symmetric with respect to the T\j = line. This symmetrization and the fact that the total polarization 
in the latter case is also smaller is due to the fact that the strength of the field is everywhere larger. In the saturated 
regime, the diagrams should be symmetric about the Ty = line, as explained above. Note also, that in the present 
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case, the polarization in the saturated regime depend on the parameter q. 

All the above considerations can be extended to more general distributions of fields (octupolar, hexapolar, etc.), 
as long as the additional components are restricted to be aligned with the dipole and hence, the global field remains 
poloidal. For a general arbitrary orientation of the quadrupole (or larger multipoles) with respect to the dipole, a 
toroidal component of the field appears, and Equation (J5J) from which all our analysis derives is no longer applicable. 
A further generalization to include these more general geometries is, nonetheless, possible along similar lines to the 
method followed here. 

7. CONCLUSIONS 

We have calculated the polarized emission by an unresolved stellar dipole, and Equations (fl7)) - (|T8| are the main 
result of this paper. To obtain them, we have derived a general expression (Equation [5]) and then made several 
approximations . 

The emissivity in a spectral line in a given direction is simply given from the atomic polarization of its upper level 
expressed in that ray reference system (Equations dU)— d3J) ) - The statistical equilibrium (balance between radiative 
processes, collisions, and magnetic precession) of the atomic polarization in the magnetic field reference system is 
given by Equation (|4]). The symmetries of the radiation field are naturally expressed by the radiation field tensors 
Jq in the atmosphere's local vertical reference system. By rotating from the local vertical to the local magnetic 
field system, then to the common dipole reference system, and finally to the LOS reference system, we get the exact 
Equation [[5]). giving the statistical tensors in a common reference system and which is then amenable to integration 
over the stellar disk for a poloidal distribution of magnetic fields. 

We have made the following approximations. 1) Assuming that the radiation field is axially symmetric with respect 
to the radial direction, we can express the disk integrated atomic polarization (hence, scattering polarization) in terms 
of just the three real functions fp(i,Td), gp(i,Fd), and hp(i,Td), whose averages over the stellar disk are shown in 
Figure 2. 2) The mean intensity Jg and the anisotropy Jq may be calculated from the intensity in a given model 
atmosphere. Alternatively, here, we estimate their value at the stellar surface (t = 1) from the emergent center-to- limb 
variation. This is the strongest approximation in our approach. The actual values of the radiation field tensors can 
be derived from more detailed theoretical calculation or observations, or be left as a convenient parametrization to be 
constrained from observations (together with the magnetic field strength and orientation). 

This approach allows us to derive explicit expressions for the distribution of intensity and polarization over a resolved 
disk and the behavior of the polarized flux for unresolved stellar dipoles. The results can be applied to unresolved 
stellar dipoles, as well as to planetary dipoles (resolved or not) at long wavelengths where the emission of the planet 
dominates over the stellar illumination. 
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